Fractional quantum Hall edge: Effect of nonlinear dispersion and edge roton 
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According to Wen's theory, a universal behavior of the fractional quantum Hall edge is expected 
at sufficiently low energies, where the dispersion of the elementary edge excitation is linear. A 
microscopic calculation shows that the actual dispersion is indeed linear at low energies, but deviates 
from linearity beyond certain energy, and also exhibits an "edge roton minimum." We determine 
the edge exponent from a microscopic approach, and find that the nonlinearity of the dispersion 
makes a surprisingly small correction to the edge exponent even at energies higher than the roton 
energy. We explain this insensitivity as arising from the fact that the energy at maximum spectral 
weight continues to show an almost linear behavior up to fairly high energies. We also formulate an 
effective field theory to describe the behavior of a reconstructed edge, taking into account multiple 
edge modes. Experimental consequences are discussed. 
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I. INTRODUCTION 

The edge of a fractional quantum Hall (FQH) sys- 
tem [1] constitutes a realization of a chiral Tomonaga- 
Luttinger liquid (CTLL). (The word chiral implies that 
all the fermions move in the same direction). In a sem- 
inal work. Wen postulated that the CTLL at the FQH 
edge is very special in that the exponent characterizing 
its long-distance, low-energy physics is a universal quan- 
tized number, which depends only on the quantized Hall 
conductance of the bulk state but not on other details 
|2|3j . He described the FQH edge through an effective 
field theory approach (EFTA) based on the postulate 
that the electron operator at the edge of the v = 1/m 
FQH state has the form 
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where (j){x) is the bosonic field operator. The imposi- 
tion of antisymmetry forces m to be an odd integer [2], 
which, in turn, leads to quantized exponents for various 
correlation functions. In particular, it predicts a relation 
/ ~ ]/3 between the current (/) and the voltage {V) for 
tunneling from a three-dimensional Fermi liquid into the 
1/3 FQH edge, which has been tested experimentally by 
Chang et al. and Grayson et al. [JHI]. 

Wen's theory describes the FQHE edge in the asymp- 
totic limit of low energies and long distances. The CTLL 
description is inapplicable at energies comparable to or 
larger than the bulk gap, where bulk excitations become 
available; we will not consider such high energies in this 
work. However, even in a range of energies below the 
bulk gap, deviations from the ideal asymptotic behav- 
ior may arise because the dispersion of the elementary 
edge excitation deviates from linearity and also exhibits 
an "edge roton minimum." The aim of this paper is to 
estimate these corrections from a microscopic approach. 
To focus on these corrections, we make appropriate ap- 
proximations (mainly a neglect of composite fermion A 



level mixing, discussed previously [8]) that guarantee an 
ideal quantized behavior at very low energies. 

A deviation from linearity in the dispersion of the el- 
ementary edge excitation is expected to produce correc- 
tions for the following reason. In the bosonic model of 
edge excitations, the spectral weights of all excitations 
at a given momentum obey a sum rule (see Eq. (|20P), 
first demonstrated by Palacios and MacDonald [9] , which 
is valid up to a unitary rotation of the basis. The long 
time behavior of the Green function and the differential 
conductance for tunneling from an external Fermi liquid 
into the FQH edge, on the other hand, are sensitive to 
the states within an energy slice. However, for linear 
dispersion, the energy and momentum are uniquely re- 
lated, so the sum rule is also valid for all states at a given 
energy, which produces a quantized power law exponent 
for the differential conductance. (A more detailed dis- 
cussion is given in Appendix B). A nonlinearity in the 
dispersion, on the other hand, produces an energy band 
for excitations, as shown, for example, in Fig. [6] below. 
In the absence of a unique relation between energy and 
momentum, the spectral weight sum rule is now valid for 
all states at a given momentum but not for all states at 
a given energy, and there is no reason to expect the same 
power law behavior as that at low energies. 

In this paper, we compute the edge spectral function 
from a microscopic approach using the method of com- 
posite fermion (CF) diagonalization, wherein we consider 
a truncated basis of states that contain no pairs of elec- 
trons with angular momentum unity. These are the only 
states that survive when the Haldane pseudopotential 
[lOlllj Vi is taken to be infinitely strong; all states con- 
taining pairs with angular momenta equal to unity are 
pushed to infinity. Laughlin's 1/3 wave function 12] is 
exact for this model. Restriction to this subspace is also 
tantamount to considering edge excitations within the 
lowest A level (or CF Landau level) . A neglect of A level 
mixing has been shown to be very accurate for the bulk 
physics, and we explicitly confirm below, by comparison 



to exact diagonalization results for small systems, that it 
provides a good first approximation for the edge excita- 
tions as well. Restricting to this truncated basis allows us 
to study systems with a large number of particles, provid- 
ing better thermodynamic estimates than were available 
previously. 

We determine the thermodynamic limit of the disper- 
sion of the elementary edge excitation. Our results show 
that while it is linear at low energies, it begins to deviate 
from linearity at an energy that is a fraction the 1/3 bulk 
gap. It also exhibits a roton minimum, which vanishes 
at a critical setback distance signaling edge reconstruc- 
tion, in agreement with previous work |13|14j . We show 
that while the spectral weights of individual excitations 
depend on various parameters |14| . they accurately obey 
the sum rule mentioned above over the entire parameter 
range that we have studied. 

To determine the effect of the nonlinear dispersion 
on the edge exponent we employ a hybrid approach de- 
scribed in Sec. V, wherein we build the spectrum from 
the elementary edge boson with a nonlinear dispersion 
but assume the spectral weights of the EFTA model. We 
evaluate the tunneling I —V characteristic and find that 
the calculated exponent remains unchanged to a very 
good approximation even at energies above the edge ro- 
ton energy where the dispersion is nonlinear. Ziilicke and 
MacDonald jT5] have also calculated the spectral func- 
tion and the I — V characteristics for a i/ = 1/3 edge 
by assuming a dispersion e(g) ~ —q In(aq) for the edge 
magnetoplasmon, where q is the momentum and a is a 
constant; they have found that the edge exponent varies 
as the inverse filling. Our calculation is based on a mag- 
netoplasmon dispersion that is obtained from a micro- 
scopic calculation for a system with Coulomb interaction 
and a realistic confinement potential. Recently, the ef- 
fect of a nonlinearity of the fermionic spectrum on the 
long-distance, low-energy correlation functions has been 
studied in Refs. |16ll7j . However, this analysis considers 
a system in which both right and left-moving modes are 
present and interact with each other, and it is not clear 
whether the same analysis would be applicable to a FQH 
edge with a single chiral mode. 

One may also expect some signature in tunnel trans- 
port that may be associated with the edge roton, which 
would then allow such transport to serve as a spectro- 
scopic probe of the edge roton. However we find that the 
effect of edge roton on tunnel transport is negligible, be- 
cause the spectral weight in the edge roton mode is very 
small. 

Our paper is organized as follows. Section H contains 
a description of the model and the method of calculation. 
In Sec. Ill, we evaluate the energy spectra for small sys- 
tems and compare them with the exact results. In Sec. 
IV, we study large systems and extract the thermody- 
namic edge dispersion, edge reconstruction and the edge 
roton. In Sec. V, we calculate the spectral weights and 
the associated sum rules for the EFTA, and we test the 
validity of these rules for the electronic spectra. We also 



outline our hybrid approach to calculate spectral func- 
tion and tunneling density of states. In Sec. VI, we 
calculate the spectral function and tunneling density of 
states, present our main results on the I — V characteris- 
tics, and mention their implications for the robustness of 
the edge exponent under a nonlinear dispersion. In Sec. 
VII, we discuss a system with a reconstructed edge using 
a field theoretic approach and address the effect on the 
tunneling exponent. We conclude in Sec. VIII with a 
summary and a discussion of the causes and implications 
of our main results. In the Appendix we give a math- 
ematical formalism for the spectral weights, sum rules, 
and the Green's functions for an ideal and a non-ideal 
EFTA. 



II. MODEL AND METHOD OF CALCULATION 

A. Hamiltonian 

We consider a two-dimensional electron system in a 
plane. The confinement is produced by a neutralizing 
background with uniformly distributed positive charge 
in a disk (denoted ^In) of radius Rn — {\/'^N/v) I; here 
N is the number of electrons, v is the filling factor, and 
I — \Jhc/eB is the magnetic length. (The symbol I is also 
used for angular momentum later, but the meaning ought 
to be clear from the context.) The background charge 
disk is separated from the electron disk by a setback dis- 
tance d. The ground state of the electron is determined 
by a microscopic calculation; we expect the electrons to 
be approximately confined to a disk of radius Rj^ to en- 
sure charge neutrality in the interior. This system is 
modeled by the following Hamiltonian: 
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where the terms on the right hand side represent 
the kinetic, electron-electron, electron-background, and 
background-background energies, respectively. Here mi, 
is the band mass of the electrons, p.- is the momentum 
operator of the jth electron and rj is its position, Aj is 
the vector potential at Tj, po = i^/27rF is the positive 
charge density spread over a disk of radius i?jv j and e is 
the dielectric constant of the background semiconductor 
material. At large magnetic fields, only the lowest Lan- 
dau level states are occupied; hence the kinetic energy 
/iaJc/2 (where uJc = eB/mi,c is the cyclotron frequency) 
is a constant which will not be considered explicitly. 



B. Electron States 

The single particle states in the n*'* Landau level are 
given, in the symmetric gauge, by 
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(3) 
where L^{x) is the associated Laguerre polynomial ^TE\ . 
n and m denote the Landau level index and angular mo- 
mentum index respectively, z = x — iy represents the 
electron coordinates in the complex plane, r = |z|, and 
all lengths are quoted in units of the magnetic length I. 
The lowest Landau level states (n = 0) are of special 
importance for our calculations below; they are given by 
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The many-body states are formed by taking lin- 
ear combinations of antisymmetric products of single 
particle wave functions denoted by |pi,p2, ■ ■ 'Pw) = 
'^Pi^L ■ ■ '"pnI*^)' "^J^ere Pj = {n^.m^i) is the single par- 
ticle state index of the «*'' electron, and a' is the cor- 
responding creation operator. We will be interested in 
the edge excitations of the FQH state at j^ = 1/3 be- 
low. The ground state has total angular momentum 
Mo = 3N{N — l)/2. The angular momentum of the 
excited state, AM, will be measured relative to Mq. 



C. Models of FQH Edge 



We have calculated the edge spectra for both smooth 
edge and sharp edge, with cutoff ^o = 2. We show in Sec. 
Ill, that the low-energy branch of the sharp edge matches 
with that of the smooth edge, and hence the edge disper- 
sion is not very sensitive to this issue. The calculations 
for the spectral functions are carried out for a smooth 
edge only. We note that a sharp edge eliminates several 
higher energy states, but does not significantly affect the 
low-energy branch and hence edge reconstruction. 



D. Exact Diagonalization 

The exact interaction energy for FQH systems can be 
calculated for small systems using numerical diagonal- 
ization techniques. In the disk geometry with symmetric 
gauge, for a given total angular momentum M, the ba- 
sis states |mi, 7712, • ■ ■ "t-jv) in the lowest Landau level are 
generated according to the conditions: 



y rrij — M: < mi < 1112 ■ ■ ■ < rriN. 



(6) 



Restricting to the lowest Landau level, the Hamiltonian 
in the second quantized representation is 



i/ = - ^ {r,s\Vcc\t,u)alalatay 



r,s,t,u 



+ '^{m\Vch\m)al^ajn + Vhh- 



(7) 



The tunneling of electrons from a two-dimensional elec- 
tron gas into a Fermi liquid (such as a metal or n+ 
doped GaAs) has been studied experimentally in two 
geometries: point-contact geometry and cleaved-edge- 
overgrowth geometry [7] . These are believed to represent 
realizations of smooth and sharp edges, respectively [TO] . 

In the point-contact geometry, the boundary of the 
two-dimensional electron gas is smooth. Theoretically 
a smooth edge can be modeled by including all possi- 
ble many-body edge excitation states for a given total 
angular momentum M {— X]i=i "^i)j placing no restric- 
tions on the maximum single particle angular momentum 
rui . The smoothness is ensured by states extending a few 
magnetic lengths beyond the disk edge. 

The cleaved-edge geometry is characterized by a long 
and thin tunneling barrier with a typical barrier width of 
about one magnetic length. Recent experiments suggest 
that the cleaved-edge-overgrowth represents the realiza- 
tion of a sharp quantum Hall edge [19] . A sharp edge can 
be modeled [13] by excluding the single particle angular 
momenta beyond a cutoff nimax, given by 



rrimax = 3(iV- 1) -I-/0, 
where Iq is taken to be a small integer. 



(5) 



Here the electron-electron and electron-background in- 
teraction matrix elements are defined as 
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with r + s = t + u. 

The background-background interaction energy per 
particle is calculated analytically to be 
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(9) 
with the energy measured in units of e^ /el. This adds 
a constant term to the matrix elements of the Hamilto- 
nian. (The constant background-background interaction 
must be included to obtain a sensible thermodynamic 
limit for the energy, but is irrelevant for energy differ- 
ences). Computing the electron-background interaction 
(K^b) requires numerical integration. To this end, we can 
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FIG. 1: (Color online) The function Fh{ri/RN\d) that ap- 
pears in the expression for the electron-background energy 
(Eq. (lOl) for several values of d and A^. Here r is the dis- 



tance of the electron from the center of the disk, and Rn is 
the radius of the disk of the neutralizing positive background 
charge. 



write the electron-background energy as 
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For d = 0, the integral in Eq. ( |10[ ) on the right hand 
side can be calculated analytically, the result of which 
has been given by Ciftja and Wexler [20]. For d ^ 0, 
numerical integration is necessary. Figure [T] shows plots 
of the function Ff, for different d and N. 



For the electron-electron interaction, we find the an- 
alytical expressions for {r, s\Vcc\t,u) given by Tsiper in 
Ref. [21J to be useful. It is then straightforward to con- 
struct the Hamiltonian matrix and diagonalize it either 
by using standard diagonalization procedures for small 
systems {N < 7) to get the full spectrum, or by the 
Lanczos algorithm for slightly larger systems {N — 8, 9) 
to get the low-energy spectrum. In the present work, we 
have performed full diagonalization for systems with up 
to 7 particles. 



E. CF diagonalization 

We exploit the fact that the CF theory produces very 
accurate wave functions for low-energy eigenstates of the 
problem. Our approach will be to construct a truncated 



basis for the wave functions for the edge excitations [11] 
[22] , and then diagonalize the full Hamiltonian within this 
basis to obtain various quantities of interest. The method 
has been described in detail in the literature [8I23I24| . so 
we present only a brief outline here. 

For the fraction u = n/{2np + 1), the CF theory 
maps interacting electrons at total angular momentum 
M to non-interacting composite fermions at M* = M — 
pN{N — 1) |25l26j by attaching 2p flux quantum to each 
electron. The ansatz wave functions **'^ for interact- 
ing electrons with angular momentum M are expressed 
in terms of the known wave functions of non-interacting 
electrons <i>*^* at M* as follows: 



j<k 



M* 



(11) 



Here a — 1,2, ••• ,D* labels the different states, T'lll 
denotes projection into the lowest LL, and D* is the di- 
mension of the CF basis. We choose p = 1 as appropriate 
for 1/ = 1/3, and restrict $^^ to states with the lowest 
kinetic energy at M* . No lowest Landau level projection 
is required for these states, as they are already in the 
lowest Landau level. 

The Landau levels at M* transform into Landau- 
like effective kinetic energy levels of composite fermions, 
called A levels. The restriction to the lowest Landau level 
at M* is equivalent to restricting composite fermions to 
their lowest A level. More accurate spectra can be ob- 
tained by allowing A level mixing and performing CF 
diagonalization (CFD) in a larger space, but that will 
not be pursued here. As will be seen below, the lowest A 
level results are sufficiently accurate for our purposes. 

The advantage of CF diagonalization is that the dimen- 
sion D* of the CF basis is much smaller than the dimen- 
sion of the full lowest Landau level Hilbert space at M; 
this allows a study of much larger systems. Table I com- 
pares the dimensions of the full Hilbert space (D) and the 
truncated CF space (D*) for 6 to 12 particles for several 
values of AM. The dimension D increases exponentially, 
approximately as I? = 10^^ exp(2A^) for large N. This 
gives Z? w 2 X 10^® for A^ = 36 particles, in dramatic con- 
trast to D* - 10-100 for < AM < 10. Of course, the 
Hilbert space reduction comes with a cost: the CF ba- 
sis functions are much more complicated than the usual 
single Slater determinant basis functions, and diagonal- 
ization of the Hamiltonian in this basis requires many 
non-trivial steps and extensive Monte Carlo. Nonethe- 
less, CF diagonalization can be, and has been, performed 
for many non-trivial cases of interest. 

We need to evaluate the matrix elements of the 
Hamiltonian in our CF basis. If ^^^(zi, Z2, • • • ,zjv) 
and \l/^(zi,Z2, ■ • • ,zn) denote two CF states at angu- 
lar momentum M, then the electron-background and 
electron-electron energy matrix elements are given by 
(^f iK^bl^f ) and (*f iVool^f )• Their evaluation re- 
quires evaluating multi-dimensional integrals, which can 
be effectively accomplished by Monte Carlo techniques 
described in the next section. The CF basis functions 
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8,033 


55,974 


403,016 


2,977,866 


22,464,381 
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8,946 


61,575 
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184,030,746 
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1,540 


9,953 


67,696 


478,025 


3,476,314 


25,879,361 


196,384,297 
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1,729 


11,044 


74,280 


519,880 


3,752,096 


27,755,663 


209,483,911 
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1,945 


12,241 


81,457 


564,945 


4,047,402 


29,753,578 


223,373,383 
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5 


2,172 


13,534 


89,162 


613,331 


4,362,833 


31,879,397 


238,091,562 
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6 


2,432 


14,950 


97,539 


665,355 


4,700,201 


34,141,000 


253,686,437 
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7 


2,702 


16,475 


106,522 


721,125 


5,060,174 


36,545,347 


270,200,645 


15 


8 


3,009 


18,138 


116,263 


780,997 


5,444,732 


39,101,065 


287,686,698 
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TABLE I: Dimension D of the full Hilbert space (used in exact diagonalization) for several values of A'^ and AM. Here AM is 
the angular momentuni measured relative to the angular momentum Mo = SN{N — l)/2 of the ground state of the u = 1/3 
FQH state. The last column gives D* , the dimension of the CF basis in the lowest A level, used in our CF diagonalization. 
The values of D* are given for sufficiently large N where they are N independent; for small N, D* may be smaller than the 
given value. 



are in general not orthogonal to each other. They can be 
orthogonahzed by the Gram- Schmidt procedure adapted 
for CF states to produce the energy spectrum as de- 
scribed in the hterature [8^,11 23 24 . Essentially, given 
the interaction matrix Vct^js = {^a\V\'^p) and the over- 
lap matrix Oa.p = {^al'^p), the energies and eigenvalues 
are obtained by diagonalizing the matrix 0~^V . 



F. Monte Carlo Methods 

Multi-dimensional integrals can be evaluated most 
effectively by the Metropolis-Hastings Monte Carlo 
(MHMC) algorithm p7H29] . For a discussion of the ap- 
plication of MHMC algorithm to quantum many-body 
systems, in particular to quantum Hall systems, we refer 
the reader to Refs. [ll|20j . For our energy calculations, 
we find it sufficient to thermalize for 100,000 iterations 
and then average over about 10 — 20 million iterations 
for each angular momentuni. For spectral weights cal- 
culations, about 200 million iterations are required for 
the eigenvector. These numbers do not vary significantly 
with N in the range of our study {N < 45) , but the com- 
putation time increases exponentially with N and AM, 
limiting our study to systems up to A^ = 45 particles for 
the energy spectrum, and N = 27 for spectral weights. 
The energies were calculated for AM = 1 — 8 and the 
spectral weights for AM = 1—4. 



III. SMALL SYSTEM STUDIES 

The CFD approach has been well tested in the past for 
the bulk physics at various filling factors and has been 
shown to capture the behavior of FQH systems accu- 
rately. Before proceeding to larger systems, we first test 
the validity of the CFD approach for the edge excitations. 

Using the CF diagonalization procedure outlined ear- 
lier, we compute the edge excitation spectra for the 
v = 1/S state for several parameters in the range d = 0- 
2.51 and AM = 0-8. Figures [2] and [3] show comparisons of 
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FIG. 2; (Color online) Comparison of CFD spectra with the 
exact spectra for the smooth edge for 6 and 7 particles, with 
and without a positive background (upper panels and lower 
panels, respectively). For the upper panels, we have taken 
d = 0. Blue triangles are the energies obtained by CF diago- 
nalization, and the adjacent '-I-' symbols (shifted along the x 
axis for clarity) are the exact energies. The high energy parts 
of the exact spectrum are not shown. All energies are quoted 
in units of e^ /el, and are measured relative to the energy of 
the ground state at AM = 0. AM is the angular momentum 
of the excitation. 



the CFD spectra with the exact spectra for smooth and 
sharp edges (exact diagonalization is possible for slightly 
larger systems for a sharp edge because of the additional 
restriction on the Fock space), demonstrating that the 
CFD approach is essentially exact. For a sharp edge, we 
have chosen the value Iq — 2 (cf. Eq. (l5])). The exact 
diagonalization results in Fig. [3] are taken from Wan et. 
al 13 . Consistent with their conclusions, we find that 
edge reconstruction occurs for d greater than a critical 
separation. 

We note that for small systems the results for sharp 
and smooth edges are not very different, as shown in 
Fig. [4J both show edge reconstruction for d > 1.51. For 
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FIG. 3: (Color online) Comparison of the CFD spectra with 
the exact spectra for a sharp edge for A'^ = 9 particles at 
v = 1/3 with d — 0.5 — 2.0 I. Blue dots indicate the en- 
ergies obtained by CF diagonalization, whereas the adjacent 
' + ' symbols (shifted along the x axis for clarity) are energies 
extracted from the figures of Ref . [T3] • 
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larger systems, as seen below, edge reconstruction occurs 
for larger d for the sharp edge as expected. 



IV. SPECTRA AND EDGE DISPERSION 



FIG. 4; (Color online) Energy spectra for smooth and sharp 
edge excitations oi a. v = 1/3 system with A*' = 9 and d = 
0.0 — 2.5 I. Blue dots are for the sharp edge, whereas the 
adjacent black diamonds (shifted along the x axis for clarity) 
are for the smooth edge. 



Having ascertained the validity of our approach from 
comparisons to exact results, we now proceed to investi- 
gate the physics in the thermodynamic limit. We study 
the edge spectra for different sizes and approach the ther- 
modynamic limit by identifying a scaling relation be- 
tween the physical momentum Sk (we take h — 1) and 
the edge angular momentum AM. The momentum is 
related to the size of the system by fc ^ r/P, where r 
is the radius of the orb ital wave function. For edge elec- 
trons, r ^ V2MI = y/2{3{N - 1) + AM)l for a system 
with i> = 1/3. This gives the momentum of the edge 
excitation to be 



Sk^ 



AM 



VGiN~l)l' 



(12) 



Henceforth we will denote (5fc, the physical momentum of 
the edge excitation, as simply k. 

Based on our edge spectra results in Fig. [5j with the 
parameters iV = 6 — 36, d = 0-2.51, and AM = 0-8, we 
make the following observations. 

i. Data Collapse: The energy spectra for differ- 
ent system sizes collapse, indicating proper scaling 
to the thermodynamic limit. The lowest branch 
in each of the four panels corresponds to the dis- 
persion of the single edge boson, for various set- 
back distances in the range d = — 2.5/. The 
data collapse to a single curve is apparent even for 



the second lowest branch, beyond which energies 
form a continuum. A few points for iV = 36 de- 
viate slightly from the common trend in the low- 
est branch, which we believe is due to convergence 
problems for larger systems. 

ii. Edge Reconstruction: For d > dc, we observe 
edge reconstruction due to competing electron- 
background energy and electron-electron interac- 
tion energy. 

iii. Nonlinearity and Edge Rotons: The lowest 
branch, though linear at low fc, eventually devi- 
ates from linearity for all d. We extract in de- 
tail the dispersion of the single boson excitation 
for various d values in Fig. [61 with polynomial fits 
shown on the plots themselves. We observe that 
the edge dispersion is nonlinear and the "linearity 
breakdown" (defined as the point at which the de- 
viation is ^20% from linear) occurs at energies in 
the range of 0.02 — 0.04(e^/e/); in experiments, this 
corresponds to the range 0.2meV to 0.4meV. For 
d < dc ~ 1.5/, the dispersion also shows a roton 
structure with the minima around fcp ~ 1.026/^^. 
The roton gap A^j is approximately 0.056(e^/e/) for 
zero setback distance, but depends on the setback 
distance and collapses at approximately dc — 1.51. 
The analytical fits for the dispersion relations are 
given in Fig. [5] 



The nonlinear dispersion and the existence of the edge 
roton he outside the assumptions of the EFTA model. In 
the next two sections we explore their effect on the edge 
exponent that is relevant to tunneling into the edge. 



V. BOSONIZATION OF FQH EDGE 

The bosonic EFTA model is based on the idea that 
the edge excitations can be mapped into excitations of a 
bosonic system, given by 



°° b?' 



i{n,})=n^io) 



1=0 



ni\ 



(13) 



where n; is the number of bosons in the orbital with 
angular momentum I. For a given state {ni\, the total 
angular momentum and total energy are given by 



AA/ 



^Ini, 



E 



niei. 



(14) 



Furthermore, the electron field operator at filling factor 
u = l/m is given by 



^t(0) 



~i\/rn(l}(x) 



/Tje 



-i-ym(p-i-{0) —i-ymtf)^ (9) 



, (15) 



where y^ is a normalization factor. The fields (j)+{0) and 
<f>^ (0) can be expanded in terms of bosonic creation and 
annihilation operators bi and bj as 



w 



i>o ^'■ 

^-W =- T.^M~'''- (16) 



A. Electronic and Bosonic edge spectra 

We first ask if the excitation spectrum of the electronic 
problem conforms to the bosonic prediction, in which all 
excitations are created from a single branch of bosons. 
Following Ref. [O], we identify the lowest energy state 
at each angular momentum AM in the electronic spec- 
trum with a single boson excitation at I = AM, i.e., 
ni = Si^/\M- This gives the energy dispersion e; of the 
single boson state as a function of ^, where we measure 
the energy ei with respect to the energy at AM = 
(M = Mq). Using the equations 'Y^ilni = AM and 
E{ni} — X); "■'^^ tlis energies of all the bosonic states 
{n/} can be obtained and identified with the energies of 
the corresponding electronic states. We note that in our 
truncated basis, the numbers of CF and bosonic states 
are equal at each AM. 



In Fig. [TJ we compare the bosonic excitation spec- 
trum obtained in this manner with the electronic spec- 
tra computed through CFD for the edges for the cases 
iV = 9,45 and d = 0.0. The CFD spectra are shown 
in blue circles and the bosonic spectra are shown in red 
triangles. In all cases, the spectra obtained from the 
bosonic picture, with the single boson dispersion as an 
input, show a close resemblance to the electronic spec- 
tra, confirming the bosonic picture as well as the inter- 
pretation of the lowest branch as the single boson branch. 
(The bosonic description becomes less accurate with in- 
creasing N or AM, but still remains accurate for the 
low-energy states). 



B. Spectral Weights 

The relation between the electron and the boson oper- 
ators given in Wen's ansatz in Eq. (fll) leads to a precise 
prediction for the matrix elements of the electron field 
operator. We will study, following Palacios and Mac- 
Donald [Ref. |9], these matrix elements, called spectral 
weights, defined by 



<^{"a 






m 



where |{«i}) represents the bosonic state with occupation 
{n/}, |0) is the vacuum state with zero bosons, V'^(^) 
is the electron creation operator at position (with one 
dimension wrapped into a circle) , and I denotes the single 
boson angular momentum. 

Using Eqs. ([I]), ( [l6| and ([l3|, it is straightforward to 
obtain the EFTA predictions for the spectral weights 



IC 



m 



"l+n2+" 



{ni}\ 



ni!n2!---l"i2"2 



(18) 



We note that the denominator in Eq. (17) eliminates the 
unknown normalization constant y^ in Eq. (|15|) 



To obtain the spectral weights from our electronic 
spectra, we need to identify a "dictionary" between the 
bosonic states and the electronic states. It is natural to 
identify the vacuum state |0) with the ground state of 
interacting electrons aX v = l/m, denoted by |^^). The 
field operator has the standard meaning of 



^n^) = E'?'*(^W^Ev'/(^)' 



(19) 



where a[ and a; are the creation and annihilation opera- 
tors for an electron in the angular momentum I state, the 
wave function for which is given in Eq. Q. The wave 
function '^^^}{{zi\) is the electronic counterpart of the 
bosonic state |{n;}) obtained through CF diagonaliza- 
tion. Using these definitions we calculate the electronic 
spectral weights. The details of the mapping and cal- 
culational method have been discussed in a previously 
published work [2]. 
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FIG. 5: (Color online) Edge spectra as a function of the physical momentum (see Eq. (12l) for A*' = 6 — 36 particles. The 
setback distance d = 0.0 — 2.51. Data collapse for the lowest spectral branch can be seen in all the panels. Lower panels with 
d> 1.51 show edge reconstruction. 



C. Spectral Weight Sum rules 



D. A Hybrid Model 



As seen below, a sum rule for the spectral weights plays 
an important role. For ly ~ 1/m, in the bosonic EFTA, 
the sum of the squared spectral weights (SSW) is given 
by (see Appendix A for a derivation), 



sswiir 



{AM - 



-D! 



J^ini = 



AM!(m-l)! ' 
AM. 



(20) 



It is natural to ask whether the above relation holds for 
the real FQH edge. We test the validity of the sum rules 
for v = 1/3 in our model of a FQH edge by computing the 
spectral weights for system sizes iV = 9 — 27 and AM = 
1 — 3. The results for individual spectral weights have 
been published in a previous work by two of the authors 
[14] . In Fig. |8J we show the plots of the SSW for Coulomb 
interactions for different N. The thermodynamic limit 
for the SSW approaches the expected result according to 
Eq. (pOl). 



To obtain results for the spectral function and the tun- 
neling density of states in the parameter regime of our 
interest, bigger systems and larger angular momenta are 
needed. We have found that it is computationally infea- 
sible to calculate the spectra for TV > 50 and AM > 8, 
and spectral weights for TV > 27 and AM > 4. To make 
further progress we used a hybrid approach. We work 
with the single boson dispersion obtained by the micro- 
scopic theory, but we assume that (i) the full spectrum 
can be constructed from it by assuming that the bosons 
are noninteracting, and (ii) the spectral weights of indi- 
vidual states are given by the EFTA model. With these 
assumptions, our model tests only the effect of nonlin- 
earity of the single boson dispersion. Corrections to the 
edge exponent arising from coupling to states outside of 
our restricted basis, as well as those from a redistribu- 
tion of the spectral weights between states, are outside 
the scope of our present study. 

As an illustration of our hybrid approach, we have plot- 
ted in panel (c) of Fig. [9] the spectral weights of various 
excited states discussed in Sec. |VA| The figure iUustrates 
that the spectral weights corresponding to a given num- 
ber of bosons have roughly the same energy, in agreement 
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FIG. 6: (Color online) Dispersion e{k) of single edge boson for different setback distances in the range d = — 1.51. The solid 
lines are a fifth order polynomial fit to the lowest branch of the energy spectra in Fig. |5] Arrows indicate the energy beyond 
which the dispersion becomes nonlinear, as defined in section IV. Dispersion minima for the un-reconstructed edges {d < 1.51) 
have been fitted with a roton curve eroton{k) = A + b{k — fco)^. Roton gap, momenta and curvature corresponding to the roton 
minima are also shown. Note that Fig. [5] has d up to 2.51, but in this figure we show the plots for the relevant distance range 
d = 0- 1.51. 
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FIG. 7: (Color online) Energy spectrum for the edge excita- 
tions of z/ = 1/3 (blue dots) , obtained by CF diagonalization, 
for A'" = 9, 45 at d = 0.0. Red triangles (shifted along the x 
axis for clarity) show the bosonic spectra generated from its 
lowest branch (see section V.A for explanation). 



VI. SPECTRAL FUNCTION AND TUNNELING 
DENSITY OF STATES 



The positive energy part of the electron spectral func- 
tion is given by [ll|30j . 

A>(A;,£;)^^|(a,iV + l|4.|0,7V)|^<5(i?-i?r'+O, 

a 

(21) 
where a denotes many-body energy eigenstates, and k, c|. 
denote the momentum (or any other) quantum number 
and the corresponding electron creation operator respec- 
tively. For the FQH edge, if we restrict to the states in 
the lowest Landau level, '^^^}{{zi}) would correspond 
to \a,N + 1). Using the definition of the spectral weight, 
we can write the spectral function as 



with previous work by Ziilicke and MacDonald 
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y^lni = A^ fc. 



(22) 
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FIG. 8: (Color online) Sum of spectral weights (Eq. ( |18[ )) obtained from the electronic spectra (see section V.B for definition) 
at angular momenta AM = 1 — 3 and d = — 2.0/. For A A/ — 1, there is only one one state, which is independent of d in our 
model. For other cases we find that the spectral weight sum is independent of d. Also, the thermodynamic limit is consistent 
with the sum rule derived from the EFTA (see Eq. po])). 



where |C{„,}p is the electronic spectral weight, E^_^^} 
is the energy of the electronic spectra measured from 
the ground state of A^ + 1 particles, and e is the en- 
ergy of the edge excitation measured with respect to the 
chemical potential /i. He re, 77 is the normalization factor 
|(0|^1'(6')|0)p in Eq. ^. We divide the energy into dis- 
crete bins [e — 6/2, e + 6/2) of width 6 and sum over the 
spectral weights for states with the corresponding ener- 
gies and momentum k to calculate A(k, e). As discussed 
in Sec VD[ we have used the electronic energy dispersion 
and the bosonic spectral weights to calculate the spectral 
function. 

In Fig. p\ panel (a), we show the energy spectra with 
spectral weights (colored) for bosonic states. The low- 
energy states have comparatively smaller weight. The 
spectral function A{k, e) (unnormalized and in arbitrary 
units) for different momenta k are shown in panel (b). In 
panel (d) we note that the energy corresponding to the 
maximum of spectral function closely follows the line of 
maximum energy for a given momenta. 



A. Tunneling density of states 

When an electron tunnels between two weakly cou- 
pled systems (labeled L, R) with a chemical potential 
difference eV , the tunneling current can be shown to be 
(PEQ]), 



I{eV) 



a,0 



reV 

r„,0|2 / dEALia,E)AR{l3,eV-E), 
Jo 



(23) 

where a, /3 are the quantum numbers of the electron 
states in the two systems, and T^.^ is the matrix ele- 
ment connecting the two states. If the energy range of 
tunneling is small, Ta^/3 can be approximated by a con- 
stant T independent of the quantum numbers. Further 
assuming one system (say L) is a metal, whose the den- 
sity of states is almost constant near the Fermi surface, 
gives the differential conductance as proportional to the 



tunneling density of states in the other system (i?, la- 
beled as "edge") as 



dV 



metal— cdsc 



For a FQH edge, using Eq. (22 1, the tunneling density 
of states (the superscript A^ -t- 1 is omitted for brevity) is 
given by, 

i?edgc(e) ^ E \CM\'Si^-EM)- (25) 

{n,} 

The relation between / and V is given by 

lieV) - J2 r deA>{k,e) 
k -^0 

~ r A^|C{„jpJ(e-i?{„,}), (26) 

which is essentially the sum over all the squared spectral 
weights of states with excitation energy e < eV (Ref. 

In Fig. 10 we show the I—V characteristics computed 
for a system of A^ = 75 particles. Log-log plots in these 
panels show several plateaus and steps in the low voltage 
region, which are purely due to the finite size effect of 
summing over a discrete set of spectral weights (in the 
low-energy regime we have very few states in spite of the 
fairly large number of particles considered). We observe, 
surprisingly in view of the physics described in the in- 
troduction, that the exponent a in / ~^ V^ remains very 
close to the ideal EFTA result of 3, within numerical er- 
rors. To explore the reasons behind the robustness of the 
edge exponent to nonlinearities in the dispersion, we have 
plotted the energy at the maxima of the spectral func- 
tion A{k, e) as a function of k in panel (d) of Fig. l9] In 
the energy region of interest (e < 0.2), the peaks roughly 
follow the ideal EFTA line. The low-energy states near 
the lower edge of the dispersion have comparatively less 
spectral weight, and their contribution to the tunneling 
density of states is negligible. 
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FIG. 9: (Color online) (a) Bosonic spectra generated through electronic dispersion for TV = 36 particles and d — 0.0 using 
the hybrid approach described in Sec. V.D; the bosonic spectral weights are shown in graduated colors, (b) Spectral function 
calculated from Eq. ( |22[ | for different k using the data in panel (a) . (c) Bosonic spectral weights plotted as a function of the 
energy, grouped by the number of bosons, (d) Energy at the maxima of the spectral function for different k. The black empty 
squares indicate the situation for a linear dispersion corresponding to EFTA. The system has A'^ = 36 particles; the setback 
distance is d = 0.0; and we have restricted to angular momentum up to AMmax ~ 19. 



B. Irrelevance of edge roton in tunneling 



VII. EFFECTIVE APPROACH FOR 

RECONSTRUCTED EDGE AND TUNNELING 

EXPONENT 



One might ask whether the edge roton produces any 
signature in a tunneling experiment. In panel (a) of Fig. 
|10[ no significant structure is seen when eV is equal to the 
roton energy. Panel (b) corresponds to the setback dis- 
tance where the roton gap just vanishes. Again, there is 
no prominently visible structure that may be attributed 
to the roton energy. An increase in the density of states 
at very low energies, shown in the log-log plots, is associ- 
ated with the edge roton, but such a signature would be 
difficult to detect in experiments. We surmise that the 
spectral weight in the roton mode is too small for it to 
be observable in tunnel transport. 



For systems which undergo edge reconstruction (in our 
case d > dc), a logical procedure would be to study the 
excitations around the new ground state which now oc- 
curs at a finite AM. This, however, is not possible in our 
numerical calculations because of computational limita- 
tions. For example, in the 27 or 45 particle system, we 
cannot go to large enough values of AM to identify the 
minimum energy. 

To make further progress, we make the assumption 
that the edge-reconstructed system can be described by 
multiple chiral edges (we take three chiral edges below), 
which interact with one another. For want of a better 
description, we further assume that each chiral edge can 
be modeled by the EFTA Lagrangian and ask to what 
extent this can describe the experiments. 

We will use the technique of bosonization to study the 
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FIG. 10: (Color online) (a) The inset shows / as a function of V for tunnel transport between a FQH edge and a Fermi liquid 
computed from Eq. ( 26 1 using the numerical edge dispersion for N = 75 particles. The main panel shows the log-log plot 
oi I — V characteristics, to better bring out the power law behavior. The red line marks the curve I cc V^. The step like 
deviations at small energies are a finite size artifact (see text for explanation). A logarithmically small deviation can be seen 
at eV = Aroton, the position of which has been marked, (b) Same as in panel (a) but for a setback distance where the roton 
gap A vanishes. The resulting slope is still close to 3 at energies e < 0.15. 



effects of density-density interactions between three chi- 
ral modes. A model with three modes (two of them mov- 
ing in one direction and the other mode moving in the 
opposite direction) is motivated by the analysis in Refs. 
[31I32J . We denote the modes as 1, 2 and 3, of which 
1 and 3 move to the right (from x = — oo to oo) and 2 
moves to the left (from a; = oo to — oo). We take the 
bosonic Lagrangian density for the system to be of the 
EFTA form (we use a slightly different normalization for 
the bosonic field in this section than in Sec. V), 



tion at time i = 0, 
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47r 



3 

dx[^£pdt(l)pdx(l)p 
p=i 



y^ dx(t>pKpqdx(l 



(27) 
where Cp — \ and —1 for the right and left moving fields 
respectively, and Kpq is a real symmetric matrix whose 
off-diagonal entries give the strengths of the interactions 
between pairs of modes; we have assumed these interac- 
tions to be short-ranged for simplicity. We have absorbed 
the velocities Vp of the modes in the diagonal parameters 
Kpp. For repulsive density-density interactions, the di- 
agonal elements Kpp as well as Ki^ are positive, while 
Ki2 and K2'i are negative; this is because the densities 
of fields 1 and 3 are given by {\/2TT)dx4'i and {l/2'K)dx4>^, 
while the density of field 2 is —{\/2-K)dx4'2- Note that the 
filling factors Vp {p = 1,2,3) have not been introduced 
in the Lagrangian density in Eq. ( |27[ ). They will appear 
later when we consider the electron creation operator on 
edge p, namely, ipp ^ e 



ict>p/^ 



To quantize the above theory, we impose the equal- 
time commutation relations [(f>p{x), Pq{y)] = —iSpqS{x — 
y). These are satisfied if the fields have the decomposi- 



dk 



<f>p{x) = / -ribpkt 






+ Kk^' 






(28) 



where [bpk, 6^^,,] = 6pq k S{k - k'). 

In the absence of off-diagonal interactions (i.e., Kpq — 
if p 7^ q), the velocities of the three modes are given 
by Kii, —K22 and K^^] the first and third are positive, 
while the second one is negative. When the off-diagonal 
interactions are present, the Lagrangian density in Eq. 
(271 can be diagonalized either by a Bogoliubov trans- 
5] or, equivalently, by solving the equations 
In the latter method, we assume that the 



formation 

of motion. 

fields take the form 



X„„e*'^'^ "°*\ where the in- 



dex a (= 1, 2, 3) labels the three different solutions, and 
Va denote the corresponding velocities. The equations of 
motion then show that the eigenvectors Xpa (which are 
real) and the velocities Va are solutions of the equations 



9=1 



(29) 



We again assume that the new velocities Va are positive 
for a = 1, 3 and negative for a = 2, and we define Cq, = 1 
for p = 1, 3 and —1 for a = 2. The eigenvectors Xpa can 
be normalized so that 

3 

2_^^p^aXpaXpfi — 5ap, 
p=l 

3 
/ ^ ^p^Oi-^pa-^qa -- ^pq- V'^^J 
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If bpk and b^k denote the original and new (Bogoliubov 
transformed) bosonic annihilation operators, we find that 
these are related as 

bak = '^Xpa[-{l + epea)bpk- -{l-epea)blf^], 

3 1 ~ 1 

bpk = '^Xpa[-{l + epea)bak + 7^(l-epea)bl,^]. 



a=l 



(31) 



Using Eq. (30 1, we can verify that [bpk,bqk'] = and 



[bpk, ^gfc'] = ^pq k S{k - k') imply that [bak, bpk'] = and 



\bak, bpk'] ~ ^aji k 5{k — fc'), as desired. 

Let us now consider the electron creation operator on 
one of the three edges, say, i/'i ~ ^i-'f'i/V^ ^ where we 
have assumed that edge 1 is associated with the filling 
factor vi. In the absence of the off-diagonal interactions, 
-01 has the scaling dimension l/(2i/i). In the presence 
of interactions, we find from Eq. (31) that the scaling 
dimension of f/'i is given by 



di 



{Xi 



{X,2? 



{Xi,f 



2vi 



(32) 



Since the second equation in Eq. ( 30 ) with p — q — \ 



implies that (-'i^ii)^ — (-^^12)^ + (-'^I'sF — !> the expres- 



sion in Eq. (32) is larger than \/{2ui) if X12 7^ 0, i.e., if 



there is a non-zero interaction K12 between modes 1 and 
2. We thus see that interactions between two counter- 
propagating modes lead to an increase in the scaling di- 
mension of the electron operator; hence the exponent for 
the two-point correlation function for electrons becomes 
larger than l/i^i. In particular, if an edge corresponds to 
a filling factor of 1/3 or less, the electron correlation ex- 
ponent on that edge will be larger than 3. Thus a model 
with multiple chiral modes in which counter-propagating 
edges interact with each other has difficulty in explaining 
the results of tunneling experiments [3HT which measure 
an exponent of about 2.7. 

Wan et al |T3j and Joglekar et al 34 also studied the 
reconstruction of FQH edges at 1/ = 1/3 and showed that 
the presence of counter-propagating edges leads to a non- 
universal exponent. Yang [31j introduced an action which 
has cubic and quartic terms in bosonic fields and showed 
that this leads to an exponent slightly larger than 3. In 
contrast, we have considered a standard action which is 
quadratic in bosonic fields and have shown that inter- 
actions between counter-propagating modes necessarily 
leads to an exponent larger than 3. 



dispersion and the associated bosonic spectra, and the 
use of spectral weights from the bosonic theory. 
The conclusions of our work are as follows. 

i. The edge dispersion is linear for energies below 
0.02-0.04e2/e/ (0.2meV to 0.4meV) depending on 
the electron-background separation. For d < dc = 
1.5/, an edge magnetoroton is observed. The maxi- 
mum roton gap is A « 0.056(e^/e/) for zero setback 
distance. 

ii. Edge reconstruction occurs beyond a critical 
electron-background separation dc ~ 1.5/ for 
smooth edges of a j/ = 1/3 system, in agreement 
with the previous literature [Oj . 

iii. A bosonic description of the edge excitation spec- 
trum is satisfactory. It requires the dispersion of 
the single boson excitation as an input. 

iv. The spectral weights of the electronic dispersion, 
though individually different from that of predic- 
tions of the bosonic theory, obey the same sum 
rules for a given angular momentum (provided A 
level mixing is neglected). 

V. The tunneling exponent is surprisingly insensitive 
to the nonlinearity in the edge boson dispersion. 
The peaks of the spectral function for different 
momenta roughly follow the linearity of the ideal 
EFTA. The low-energy states have a small spectral 
weights and contribute negligibly to the tunneling. 

vi. The roton has no significant contribution to the 
spectral function and hence to the tunneling den- 
sity of states. Only a logarithmically weak signa- 
ture of the roton may be observed in tunneling ex- 
periments. 

vii. It is well known that the model assuming a single 
chiral mode is not adequate for understanding the 
results of experiments on systems which undergo 
edge reconstruction. An effective theory descrip- 
tion with three chiral edges sX v — 1/3 produces 
an exponent that is larger than 3, contrary to the 
experimental finding of a smaller-than-3 exponent. 
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X. APPENDIX 



for the ID chiral edge is 



A. Sum rules 



The derivation of the sum rules in Eq. ( pO| for squared 
spectral weights at a given angular momenta X]|n > '"' ~ 
M is given below. Consider the multinomial expansion 
(Ref. [IF p. 823) 



OO , I. 



E 

\k=l 



EE 



{aj} 






E-?'"j " "' 



(33) 



G{x,t) = (0|T(*(x,i)*^(0,0))|0) 

= (0|e-*-^*e''=°*(0,0)e-''=="e'^**^(0,0)|0), t > 0. 

(38) 
We map the edge to a disk by setting x = R9 with radius 
R = 1 and insert a complete set of states within the 
subspace of single boson modes: 



^^|Af,{n,})(M,{nz}|=/,; Af = ^Zn,. (39) 

M {n,} I 



With the following transformations 

m -^ h number of bosons, 

Xfc — >■ TO inverse filling factor, 

^3 ~^ "•j bosons occupation, 

n — )■ Af angular momentum, 



we obtain 

it 



OO J, 



\k=l 



' ^-^ ^-^ M! nin,!7"0 



= ^!EE 



E-?"j = ^^■ 



(34) 



Hence we get 



b=0 \ /c = l 



OO OO 



EEE^^n 



b=0M=6{nj} i 

We simplify this by noting the relation 



jij Ij"^ 



(35) 



exp m y^ 



— = g-mln(l-t) ^ ^ 



fc=l 



(1 - ty 



(36) 



The sum of the squared spectral weights W is the coeffi- 
cient of t*^ 



E i^i:; 

in.} 



(„) 2 _ (M + m-1)! 



M!(m-1) 



'-^, ^?n, = M. (37) 



We make the following substitutions 



k = AM; \ = M/y/6{N-l) 
*^{"i} ^ ^e;n; = Au_F^?n; = Ai;_fM, (40) 



; 



I 



and proceed to calculate the Green's function: 



Gix,t) 



e --(U|*(0,0)e^^*|A/,H}) 



x(Af,{nJ|M/t(0,0)|0) 
= ^^e-*^^("-"^*)(Af,{n/}|*^(0,0)|0)|2 

M {«,} 

= ^e-'^*^("~'"'*)^|(Af,{n,}|^t(o^O)|0)|2 

{ni} 

(x-VFt)fM + m-V 

771—1 



Af 



V^ ~i\M{x-VFt) 



M 



{x — vpty 



(41) 



This shows how the power law follows from a combina- 
tion of the linear dispersion and the sum rule. For a 
general dispersion w^, we evaluate the commutators of 
the bosonic fields to find the Green's function for the 
edge: 



B. Green's function 

The ideal EFTA assumes a linear dispersion e{k) = 
vpk, and the sum rule in Eq. (37). The Green's function 



27r 



G[X,t) = -—e^^k I'" "" . 



(42) 
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